Quasiparticle states of the Hubbard model near the Fermi level 
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The spectra of the t-U and t-t'-U Hubbard models are investigated in the one-loop approximation 
for different values of the electron filling. It is shown that the four-band structure which is inherent 
in the case of half-filling and low temperatures persists also for some excess or deficiency of electrons. 
Besides, with some departure from half-filling an additional narrow band of quasiparticle states arises 
near the Fermi level. The dispersion of the band, its bandwidth and the variation with filling are 
close to those of the spin-polaron band of the t-J model. For moderate doping spectral intensities 
in the new band and in one of the inner bands of the four-band structure decrease as the Fermi level 
is approached which leads to the appearance of a pseudogap in the spectrum. 

PACS numbers: 71.10.Fd, 71.27.+a 



I. INTRODUCTION 

In the last two decades the discovery of high-T c super- 
conductors, heavy-fermion compounds and organic con- 
ductors has revived interest in strongly correlated elec- 
tron systems. One of the simplest and still realistic mod- 
els in this field is the one-band Hubbard model two- 
dimensional version of which has been extensively studied 
in connection with the cuprate perovskite superconduc- 
tors. Along with Monte-Carlo simulations, 1 '- different 
cluster mcthodsji^ 4 * 5 . the operator projection technique, 6 
the generating functional approach^ various versions of 
the diagram technique2i2ii2iiLi2ii£ are used for the inves- 
tigation of the model. In a system in which the Coulomb 
interaction dominates, it is reasonable to treat this inter- 
action exactly and the kinetic energy in the framework 
of a perturbation theory. 

In the diagram technique proposed in Refs. Illll2ll3l 
the power expansion is expressed in terms of site cu- 
mulants of electron creation and annihilation operators. 
From the expansion the Larkin equation for the elec- 
tron Green's function can be derived^ 1 * 1 ^* 1 ^ The equa- 
tion can be solved in the one-loop approximation. How- 
ever, the obtained solution has a flaw - a negative spec- 
tral weight in two narrow frequency regions^ This flaw 
can be remedied by an interpolation using results for 
regular regions^ The obtained spectral function 1 ^ was 
shown to be in agreement with results of Monte-Carlo 
simulations 2 ^ at half-filling and for moderate tempera- 
tures when the magnetic correlation length is comparable 
with the intersite distance. In particular, it was shown 
that in agreement with results of Monte-Carlo and cluster 
methods 1 ' 2 ' 3,4,5 the diagram technique is able to describe 
the four-band structure of the spectrum at half-filling. 

In the present paper the same method is used for 
the investigation of the energy spectrum at a depar- 
ture from half-filling. It is shown that in this case the 
above-mentioned four bands persist and additionally a 
new band arises in some vicinity of the Fermi level. By 
its properties - the dispersion, bandwidth and the vari- 
ation with filling - the band resembles the spin-polaron 
band of the t-J modelii 5 - For moderate doping the spec- 



tral intensities in the new band and in one of the inner 
bands of the four-band structure decrease as the Fermi 
level is approached. This produces a pseudogap near the 
Fermi level. For the hole-doped case, n < 1, the mag- 
nitude of the pseudogap observed in photoemission de- 
creases with increasing the hole doping 1 — n, while for 
the electron-doped case, n > 1, this magnitude increases 
with increasing the electron doping n — 1. Here n is the 
electron concentration. Together with the t-U Hubbard 
model the t-t'-U model is also considered for the ratio 
t 1 ft = —0.3 of the next-nearest and nearest neighbor 
hopping constants^ As for the case of half-filling the 
calculated spectral functions and dispersions appear to 
be similar to those obtained by Monte-Carlo and cluster 
methods, provided that doping or temperature are high 
enough to ensure a short magnetic correlation length. 

Main formulas used in the calculations are given in the 
following section (the detailed derivation of these formu- 
las can be found in Ref. X\M . The discussion of the ob- 
tained results and their comparison with results of other 
methods are carried out in Sec. Ill and IV for the t-U 
and t-t'-U models, respectively. Concluding remarks are 
presented in Sec. V. 



II. MAIN FORMULAS 



The Hamiltonian of the Hubbard model reads 
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where t nm is the hopping constants, the operator a* aa 
creates an electron on the site n of the two-dimensional 
square lattice with the spin projection a = ±1, U is 
the on-site Coulomb repulsion and the electron number 
operator n na — aj, ff a nff . 

The diagram technique proposed in Refs. [Tl| and[l2 is 
used in the present work for calculating Green's function 



G(nV,nr) = (To n v(T'K,(r)>, 
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where the angular brackets denote the statistical averag- 
ing with the Hamiltonian H. = H — \x ^2 na n ncr , \i is the 
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FIG. 1: Irreducible diagrams of the first four orders of the 
expansion in powers of i n m. 



chemical potential, T is the time-ordering operator which 
arranges other operators from right to left in ascend- 
ing order of times r, a na (T) — exp(7ir)a ncr exp(— Hr), 
and a no -(r) — exp(7ir)aj 1(T exp(— TLt). With the use of 
the diagram technique the following Larkin equation is 
derivecUi^ii* for the Fourier transform of function (J2J : 
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where f k = E n e "' kKm)i nm, u>i = (21 + 1)ttT is the 
M atsubara frequency with the temperature T and an in- 
teger I, K(k,ioji) is the sum of all irreducible diagrams 

- the diagrams which cannot be divided into two parts 
by cutting a hopping line. Such diagrams which appear 
in the first four orders of the perturbation expansion in 
powers of t nm are shown in Fig. ^ with their signs and 
prefactors. Here circles are cumulants 1 - of electron oper- 
ators. In diagrams, the order of a cumulant is determined 
by a number of incoming or outgoing hopping lines (di- 
rected lines in Fig.^l. Cumulants of the first and second 
orders which will be used below read 

K\ (tV, t<t) = (Ta ff (T')o ff (T)) C', 
K 2 (T'a,T<T,T{ai,Tiai) = 

(Too-CrOoo-Mao-j (r[)a ai (n)) 

-Ki(r'a, tu)Ki (tJcti, t\cfi) 

+Ki(r'a, Tiai)Ki(r[ai,Ta), 

where the subscript "0" near the angular bracket in- 
dicates that the averaging and time dependencies of 
the operators are determined by the site Hamiltonian 
H n = X3 CT [(^/2)n n cr^n,-cr — M n ncr]- All operators in the 
cumulants belong to the same lattice site. Due to the 
translational symmetry the cumulants do not depend on 
the site index which is therefore omitted in the above 
equations. 

Partial summation is implied in the diagrams in Fig. ^ 

- the irreducible diagrams are included in the hopping 
lines which therefore correspond to the expression 
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In the one-loop approximation used below, the total 
collection of irreducible diagrams K (k, iui) is substituted 
by the sum of the two diagrams (a) and (b) in Fig. ^ 
Thus, 



K(iwi) = Ki(iwj) — TyiK2(iui(r,iwi 1 (Ti,iwi 1 (Ti) 
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where Ki(iu)i) and K2(iuJia,iuJi 1 ai,iu)i 1 ai) are the 
Fourier transforms of the cumulants of the first and sec- 
ond orders, respectively, N is the number of sites and I 
set ^ k it = 0. Notice that in this approximation K does 
not depend on momentum. The cumulants reacU* 
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-gMmJ) + goi{iuJi 1 )g 12 {iuJi 1 )(gl 1 (iuji) 

+gUi Ul ))]}-Z^U^5u^' 0{Eo+E2) 
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where (3 = T~ l , E = 0, E x = -/i, and E 2 = U - 2\x 
are the eigenenergies of the site Hamiltonian H n , Eij — 
Ei-Ej, Z Q = e -P Eo +2e- f3El +e- l3E2 is the site partition 
function, gij(iui) = (iui +E ij )- 1 . 

Equations |[BJ can be significantly simplified for the 
case of principal interest U 3> T. In this case if /i satisfies 
the condition 



A < it < U - A, 



(7) 



where A > T, the exponent e~^ El is much larger than 
e -fiE anc L e -f3E 2 ^ gy p ass i n g to real frequencies one 

can ascertain that terms in ^^ K2 with the two latter 
multipliers contain the same peculiarities as terms with 
e -/3£i. Therefore terms with e~^ E ° and e _/3£2 can be 
omitted and Eq. JBJ is simplified to 
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Further simplification can be achieved by using the 
Hubbard-I approximation^ for the Green's function 
G(k, iu)^) on the right-hand side of Eq. 0. The re- 
spective expression is derived from Eq. (J3J) if the total 
irreducible part K(k, iwi) is approximated by the first cu- 
mulant K\(iwi) [the diagram (a) in Fig.p from Eqs. © 
or Q^ 9 ' 11 ! 13 This gives 
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After carrying out the summation over l\ in Eq. (J5J) 
with the use of Eqs. (f5f) — (|1(J|) it is convenient to turn to 
real frequencies by substituting iui with z = w+irj where 
r\ is a small positive constant which affords an artificial 
broadening. 



III. THE t-U MODEL 

At first let us consider the t-U model in which only 
the nearest neighbor hopping constant t is nonzero and 
tu = 2t[cos(k x ) + cos(ky)]. Here the intersite distance 
is taken as the unit of length. Due to the electron-hole 
symmetry in this case the consideration can be restricted 
to the range of the chemical potentials /i < (7/2. 

Figure[2]demonstrates Im_K(uj) calculated with the use 
of Eqs. © and IjEfy- lfTHjI . The change to real frequencies 
converts the Matsubara function (J2J) into the retarded 
Green's function^ It is an analytic function in the upper 
half-plane which requires that ImK(u)) be negative. As 
seen from Fig. this condition is violated at ujd = —fj, 
and U — \i. The problem is connected with divergen- 
cies at these frequencies introduced by functions goi(^) 
and gn{oj) in the above formulas. As can be seen from 
the procedure of calculating the cumulants, these func- 
tions and divergencies with sign-changing residues will 
appear in all orders of the perturbation theory. It can be 
expected that in the entire series the divergencies com- 
pensate each other and the resulting lmK(u>) is nega- 
tive everywhere. However, in the considered subset of 
terms such compensation does not occur. Nevertheless, 
as seen from Fig. |21 at frequencies neighboring to u>d the 
irreducible part is regular and, if the used subset of dia- 
grams is expected to give a correct estimate of the entire 
series for these frequencies, the values of ImK(ui) near 
tod can be reconstructed using an interpolation and its 
values in the regular region. 1 Examples of such inter- 
polation are given in Fig. [21 The function K(z) has to 
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FIG. 2: The imaginary part of K[u) calculated using Eqs. © 
and ©-JTUJl for a 100x100 lattice, t = -U/8 and T = 0.001(7 
(the dashed lines), (a) /x = 0.5(7, r] = 0.0117. (b) \x = 0.1(7, 
77 = 0.02(7. The solid lines show the corrected ImK(uj) (see 
text) . 



be analytic in the upper half-plane also and therefore its 
real part can be calculated from its imaginary part using 
the Kramers-Kronig relations. I used this way with the 
interpolated ImK(ui) to avoid the influence of the diver- 
gencies on ReK(cu). However, the application of the in- 
terpolation overrates somewhat values of |Im74T(oj) which 
leads to the overestimation of the tails in the real part. 
To correct this defect the interpolated K{u>) is scaled so 
that in the far tails its real part coincides with the values 
obtained from equation (J5j). 

As seen from Fig. EJl, at half-filling, /1 = (7/2, ImK(cu) 
has two broad minima. With decreasing the chemical 
potential from this value the minima shift with respect 
to the Fermi level without a noticeable change of their 
shapes until the Fermi level enters the left minimum 
which for t = —U/8 occurs at \x w 0.17(7. As this takes 
place, two new sharp minima arise near frequencies — [i 
and U — /i on the background of the above-mentioned 
broad minima (see Fig.^). The appearance of the broad 
features in Figgis connected with the third term on the 
right-hand side of Eq. ©, while the sharp minima are 
related to the second term in this formula. Its contri- 
bution to K(lo), Eq. |J5)|. grows rapidly when the Fermi 





FIG. 3: The spectral function A(kui) of the t-U model calculated for momenta along the symmetry lines of the square Brillouin 
zone in a 40x40 lattice for i = -U/8, T = 0.001(7, r\ = 0.0217, fj, = 0.2(7 (a), and fj, = 0.1(7 (b). 



level enters the broad minimum. 
The spectral function 
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obtained from such calculated irreducible part for mo- 
menta along the symmetry lines of the square Brillouin 
zone is shown in Fig. [3] The shapes of the spectral func- 
tion in Fig. |3^ are nearly the same as at half-fillingi4 - 
as in the case of K(ui), with decreasing /z from U/2 to 
0.17C7 these curves shift with respect to the Fermi level 
without perceptible changes in their shape. In agree- 
ment with results of Monte-Carlo simulations^ and clus- 
ter methods^*^ four bands can be distinguished in these 
spectra. For parameters of Fig. [3K these bands are lo- 
cated near the frequencies — 4|i|, 0, 4|t|, and 9|i|. For 
the major part of the Brillouin zone maxima forming 
the bands arise at frequencies which satisfy the equation 
1 — t^KcK(uj) = and fall into the region of a small 
damping llnxftT (w) | [see Eq. (|ll|l ]. As seen from Fig. |2 
such regions of small damping are located between and 
on the outside of the two broad minima in Im_KT(w). This 
is the reason of the existence of the four well separated 



bands - two of them are located between the minima 
of ImK (u>), while two others are on the outside of these 
minima. Broader maxima of A(kuj) for momenta near the 
boundary of the magnetic Brillouin zone are of different 
nature - since tk is small for such momenta, the resonant 
denominator in Eq. i|ll|) does not vanish and the broad 
maxima in the spectral functions reproduce the maxima 
of — \m.K{ui) in the numerator of this formula. 

More substantial changes in A(kuj) occur for /i < 
0.17C7. As seen from Fig. |2t>, the four-band structure 
persists also in this case. In addition to this there ap- 
pear sharp dispersive features near to = — /j, and U — jjl, 
the latter feature being substantially weaker in the case 
of hole doping, fi < U/2. It is clear that these changes 
in the spectral function are connected with the changes 
in ImK(u}) shown in Fig. [21 For the hole-doping case 
the peaks near —jjl axe located in the nearest vicinity 
of the Fermi level. For the parameters of Fig. |2t> nei- 
ther these peaks nor the peaks forming the lower inner 
band cross the Fermi level - their intensities decrease as 
it is approached. As a consequence a pseudogap arises 
in the spectrum near the Fermi level. Recently spectra 
with analogous pseudogaps were also obtained by cluster 
methodsi^22i2ii2^ In these works such spectral peculiar- 
ities were identified with the pseudogap observed in the 
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FIG. 4: The dispersion of maxima of the spectral function in 
the t-U model on a 20x20 lattice for t = -(7/8, T = 0.001(7, 
■q = 0.02(7, n = 0.1(7 (a), 0.05(7 (b), and 0.01(7 (c). Here 
darker areas correspond to larger intensities. The points X, 
r, M, and M' correspond to the momenta (n, 0), (0, 0), (n, ir), 
and (n/2,n/2), respectively. 



photoemission of cuprates^ 

Figure 01 demonstrates the dispersion of maxima of the 
spectral function near the Fermi level for three values 
of the chemical potential. These maxima form the lover 
inner band above the Fermi level and the new band aris- 
ing at ^w 0.17(7 below and at the Fermi level. At this 
value of the chemical potential the width of the new band 
is approximately equal to \t\ — 2 J where J = 4t 2 /U is 
the superexchange constant of the effective Heisenberg 
model which describes magnetic excitations in the limit 
U ^> |£|&i£ The bandwidth decreases with reduction in 
the electron concentration. The maximum energies of 
the band are located near the boundary of the magnetic 
Brillouin zone. The dispersion is much larger in the direc- 
tion (tt/2, 7t/2) — (0, 0) than along the boundary of the 
magnetic Brillouin zone (n, 0) — (0, n). These proper- 
ties of the band resembles those of the spin-polaron band 
of the t-J model. This latter band is also located near 
the Fermi level, has the similar dispersion and the band- 
width, which decreases with decreasing nM^ As seen 
from Fig. 01 with decreasing /i and n the Fermi level shifts 
to the lower edge of the pseudogap and enters the new 
band at \x rs 0.07(7. In the photoemission which probes 
the part of the spectral function occupied by electrons 
this change in the energy spectrum will look like the de- 
crease with the subsequent disappearance of the pseudo- 
gap with decreasing the electron concentration. Such be- 
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FIG. 5: The intensity plot of the spectral function at the 
Fermi level in the first quadrant of the Brillouin zone. The 
darker areas correspond to larger intensities. The t-U model 
on a 40x40 lattice with the parameters t — —(7/8, T = 
0.001(7, r) = 0.02(7, and /i = 0.05(7. 



havior of the pseudogap is indeed observed in hole-doped 
cuprates^ Notice that in the case of electron doping the 
pseudogap observed in photoemission will grow with n 
until the Fermi level enters the new band (see Fig. 01. 

The intensity plot of the spectral function at the Fermi 
level is shown in Fig. [5J The plot was obtained by 
averaging the spectral function in the frequency range 
[— 0.01(7, +0.01(7]. The dark area which corresponds to 
the maximum intensity can be interpreted as the Fermi 
surface in the region of the chemical potential 0.12(7 % 
jj, < 0.25(7 where the lower inner band crosses the Fermi 
level (see Fig.|3Ji), and for fi < 0.07J7 where the crossing 
occurs with the new band (see Figs.^t and c). For /j, in 
the intermediate region maxima of both bands lose their 
intensities as the Fermi level is approached (see Figs. |3d 
and|l|i). However, also in this case the intensity plot 
is similar to that shown in Fig. [SJ With decreasing /i 
and n the Fermi surface shrinks to the center of the Bril- 
louin zone and near fi = 0.01(7 changes its shape from 
a diamond centered at (%,ir) to that centered at (0,0). 
In contrast to the results of the cluster methods22*2ii22 
in the used approximation the variation of the intensity 
along the Fermi surface is small both in the t-U and t-t'-U 
models. 

For the case of half-filling comparison with the data of 
Monte-Carlo simulations 2 shows that the spectra of the 
one-loop approximation are closer to the results obtained 
at T = 0.33|t| than to those derived for T = 0.1|t|JA The 
value T = 0.33 |i| is close to the superexchange constant 
J for t = — (7/8. For such temperatures, the correla- 
tion length of the short-range antifcrromagnetic order is 
comparable to the intersite distanced Thus, it can be 
concluded that the one-loop approximation is more ap- 
propriate for short correlation lengths. This conclusion 
is also corroborated by the shapes of the spectral func- 
tion in Fig. [3^ which, as mentioned, are close to those 
at half-filling. In this figure there are no indications that 
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FIG. 6: The spectral function A(kuj) of the t-t'-U model calculated for momenta along the symmetry lines of the square 
Brillouin zone in a 40x40 lattice for t = -(7/8, t' /t = -0.3, T = 0.001(7, T) = 0.02(7, fi = 0.1(7 (a), and /i = 0.9(7 (b). 



the unit cell is doubled which manifests itself in a repli- 
cation of some parts of the quasiparticle dispersion with 
the period (tt,tt). Such a doubling is inherent in the an- 
tiferromagnetic order with a correlation length which is 
much larger than the lattice spacing. This limitation of 
the one-loop approximation is partly compensated by the 
fact, known from experiment in cuprates 26 and from the 
t-J models that the correlation length decreases rapidly 
with departure from half-filling and becomes comparable 
to the lattice spacing already at 1 — n « 0.1 even for 
temperatures T <C J. 

The spectral functions and quasiparticle dispersions 
calculated in the one-loop approximation are close to 
those obtained in Monte-Carlo simulations and cluster 
methods, provided that doping or temperature ensure a 
short magnetic correlation length (cf. Fig. [3] with Figs. 9- 
11 in Ref.Q and Fig. 2 in Ref-HH). The most important 
differences are connected with the fact that the one-loop 
approximation overestimates the spectral intensity in the 
lower inner band near the momentum (n, n) and in the 
upper inner band near (0,0). Since the Fermi level crosses 
these bands with departure from half-filling, this leads to 
underestimating (overestimating) of the electron concen- 
tration 
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at hole (electron) doping. For example, for the parame- 
ters of Figs. 12^ and b the concentrations calculated with 
the use of Eq. 1)1 2H equal to 0.87 and 0.74, respectively. 
However, from the comparison with the results of Monte- 
Carlo calculations 2 it can be concluded that the concen- 
trations have to be approximately 0.95 and 0.9, respec- 
tively. This is the reason why the spectra in the above 
figures were labeled with the chemical potential rather 
than with the concentration. 



IV. THE t-t'-U MODEL 

Now let us consider the t-t'-U model with the ini- 
tial electron dispersion tk = 2i[cos(fc x ) + cos(k y )] + 
At' cos(k x ) cos(ky). Both for the case of hole and electron 
doping the ratio t' jt = —0.3 of the hopping constants 
for the next-nearest and nearest neighbors is accepted. 16 
For both cases the spectral function and the dispersion 
of quasiparticle peaks near the Fermi level are shown in 
Figs. inland [7| As in the case of the t-U model, at half- 
filling the spectrum of the t-t'-U model contains four well 
separated bands. Also in analogy with the former model 
a dispersive feature and a pseudogap appear near the 
Fermi level at a certain level of doping. However, in the 
case of the t-t'-U model the obvious asymmetry of the 
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FIG. 7: The dispersion of maxima of the spectral function in the t-t'-U model on a 20x20 lattice for t — —U/8, t'/t = —0.3, 
T = 0.001(7, r\ = 0.02(7, ^ = 0.1(7 (a), 0.0517 (b), 0.9(7 (c), and 0.95(7 (d). Here darker areas correspond to larger intensities. 
The points X, F, M, and M' correspond to the momenta (n, 0), (0,0), (7r,7r), and (n/2,n/2), respectively. 



hole and electron doping stands out - the new feature 
and the pseudogap are less pronounced in the case of 
electron doping. This is also apparent from the compar- 
ison of Figs. [7^, and c. Such behavior is a consequence of 
the asymmetry in the filling dependence of IthK{lj) - for 
identical offsets from fi = U/2 the sharp minima which 
are similar to those shown in Fig. [2^ are more intensive 
in the case of hole doping than for electron doping. This 
asymmetry is connected with the contribution of the sec- 
ond term on the right-hand side of Eq. J^J to ImK(uj). 

Figure[7|demonstrates the dispersions of the new band 
and of the lower (upper) inner band in the case of hole 
(electron) doping, the inner band being located above 
(below) the Fermi level. The pseudogap between these 
bands becomes more pronounced with increasing h in the 
case of electron doping (cf. Figs.0: and d). Similarly to 
the t-U model, for the hole-doped case the magnitude of 
the pseudogap observed in photocmission decreases with 
increasing the hole doping (see Figs. Up, and b), while for 
the electron-doped case this magnitude increases with 
increasing the electron doping (see Figs.0; and d). 

Also as for the t-U model, shapes of the spectral func- 
tion and quasiparticle dispersions calculated in the t-t'-U 
model in the one-loop approximation are similar to those 
obtained in the cluster methods, provided that doping or 
temperature are high enough to ensure a short magnetic 
correlation length (cf. Fig. with Fig. 2 in Ref. |2(J and 



with Fig. 5a in Ref. l21fl . Again the most important differ- 
ence between these two groups of results is a larger spec- 
tral intensity in the upper inner band near (0, 0) for the 
hole-doped case and in the lower inner band near (71", n) 
for the electron-doped case in the one-loop approxima- 
tion. 



V. CONCLUSION 

In the present work the diagram technique was used 
for the investigation of the energy spectra of the t-U and 
t-t'-U Hubbard models at a departure from half-filling. 
The one-loop approximation was applied which in the 
used diagram approach is a successive improvement to 
the Hubbard-I approximation. In agreement with results 
of Monte-Carlo simulations and cluster methods at half- 
filling the obtained spectra of the models contain four 
bands. The four-band structure persists also for some 
departure from half-filling. Additionally in these condi- 
tions a new narrow band of quasiparticle states arises 
near the Fermi level. The band energy is maximum 
near the boundary of the magnetic Brillouin zone. The 
dispersion of the band is much larger in the direction 
(7r/2, 7r/2) — (0, 0) than along the boundary of the mag- 
netic Brillouin zone (tt, 0) — (0, 7r). The width of the 
band is of the order of the superexchange constant and 



s 



decreases with increasing doping. By these properties 
the new band resembles the spin-polaron band of the t-J 
model. For moderate doping the intensities of maxima in 
the new band and in one of the inner bands of the four- 
band structure decrease as the Fermi level is approached. 
As a consequence a pseudogap arises in the spectrum 
near the Fermi level. With hole doping the magnitude 
of the pseudogap observed in photoemission decreases 
and eventually the pseudogap disappears in agreement 



with experimental observations. With electron doping 
the magnitude of the photoemission pseudogap increases. 
Shapes of the spectral function and quasiparticle disper- 
sions calculated in the one-loop approximation are sim- 
ilar to those obtained by Monte-Carlo simulations and 
cluster methods, provided that doping or temperature 
are high enough to ensure a short magnetic correlation 
length. 
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